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(N ' Abstract 

In this paper the solutions of the zero-sound dispersion equation in the random phase 
approximation (RPA) are considered. The calculation of the damped zero-sound modes 
uj s (k) (complex frequency of excitation) in the nuclear matter is presented. The method 
is based on the analytical structure of the polarization operators II(a;, k). 
The solutions of two dispersion equations with II(a;,A;) and with Re(II(u;,/c)) are com- 
pared. It is shown that in the first case we obtain one- valued smooth solutions without 
"thumb-like" forms. 

Considering the giant resonances in the nuclei as zero-sound excitations we compare the 
experimental energy and escape width of the giant dipole resonance (GDR) in the nucleus 
. A with co s (k) taken at a definite wave vector k = k^. 

q 

1 Introduction 

The investigation of the zero-sound dispersion equation was started in the Landau kinetic 
theory [TJ [2]. One of the direction of the further development of this theory consisted in 
application of RPA to the problem of the zero-sound mode propagation (2j [3] . The problem 
of the propagation and spreading of the zero-sounded excitations in the various media, at the 
different temperatures and energies, is now a subject of the active investigations [4]- [8] - 

In this paper we consider the problem of zero-sound excitations in the framework of the 
standard RPA [T], EJ, E] • We reproduce the well-known results for the excitations with the real 
frequencies (for the small wave vectors k). The goal of this note is to continue calculations to 
the larger k when the zero-sound u s (k) and the quasiparticle-quasihole lp -lh modes overlap. 
In the region of overlapping we have the complex zero-sound frequencies u s (k). 

The polarization operator H(u, k) has two cuts, which are related to the appearance of 
the real p/i-pairs. The overlapping of u s (k) and p/i-pairs means that the zero-sound solution 
(its real part Keu s (k)) is placed on the cut of the polarization operator. The imaginary part 
lmu s (k) shows the shift of the solution from the real axis. There are two possibilities: 1) u s (k) 
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is shifted to the lower part of the physical sheet below the cut and 2) u s {k) is shifted to the 
unphysical sheet neighboring with the physical one through the cut. We will try to demonstrate 
that the second possibility is realized. 

The imaginary part lmu s (k) is determined by the processes which result in the appearance 
of the cut in the polarization operator. It means that the physical reason of appearance of 
lmu s (k) is the emerging of the real p/i-pairs when the real particle leaves the nucleus. Therefore 
we can interpret this Imu s (fc) as the escape width of the zero-sound excitations. 

We compare the solutions of two dispersion equations, one with the full polarization operator 
H{oj,k) and another with the real part of polarization operator Re (U(u , k)) . It is shown 
that there is an important distinction between the solutions in the region with the nonzero 
Im (H(u, k)). The dispersion equation with the full U(u,k) has the smooth solutions u s (k) 
which can be continued far enough (till the application of the RPA is valid ). The imaginary 
part of u s (k) emerges when the overlapping of the solutions and the cut of polarization operator 
takes place and is a smooth slowly decreasing function. In contrary, solutions of the dispersion 
equation with Re (n(o;, k)) have the "thumb-like" form and disappear at a special value of the 
wave vector, at this point the imaginary part go to the infinity. 

From the very beginning of the application of the Landau Fermi liquid theory to the nuclei 
[9] it was shown that the giant resonances in nuclei can be considered as zero-sound excitations 
[T0l[TT]. There are a lot of calculations of the energy and the width of the collective states and 
the giant resonances in the nuclei in the framework of the theory of the finite fermi systems [9] . 
The relation of RPA to the Landau Fermi liquid theory in the nuclei is demonstrated in j!2j . 
The theoretical consideration of the giant resonances can be characterized by some approaches. 
The first one considers the nuclei with the nuclear mean field potentials and different forms of 
the residual interaction. It uses the wave functions and the energies of the one-particle states 
with the further diagonalization of the RPA matrix jTU [13]. The second approach is related 
to the description of nuclei with the collective variables in the framework of the quantum 
hydro dynamic model [TT| IT3]. 

In this paper we use the Fermi liquid theory (9j [12] to obtain the zero-sound dispersion 
equation in the nuclear matter. Further the application of this equation solutions to the concrete 
nuclei is made. The branch of solutions of the zero-sound dispersion equation in the nuclear 
matter, co s (k), permits us to evaluate the energies and widths of the giant resonances in nuclei. 
Using the Steinwedel- Jensen model [TH [15] we can find a value of the wave vector k = kA 
corresponding to the giant resonances in the nucleus A. In the next step, ui s {kA) can be compared 
with the experimental giant resonance energies. In this work the comparison is made for the 
giant dipole resonances. The imaginary part of u) s {k) is the result of the zero-sound decay to 
the real particle-hole pairs and can be compared with the escape width of resonances. 

In Sect. 2 the explanation of our approach is given. In Sect. 3 the dispersion equation 
solutions with the full polarization operator U(u,k) and with ~Re(U(u,k)) are compared. In 
Sect. 4 we evaluate the energies and widths of the giant dipole resonances using u s {k) and 
compare the result with the known formula fitting the experimental GDR energies. 
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2 Method of calculation 



We present a method of calculation of the zero-sound solutions u s (k) of the Eq. ([I]), which are 
real or complex-valued depending on a value of the wave vector k. We use the standard RPA to 
obtain the retarded polarization operators U(cj, k) (TJ [3]. The zero-sound dispersion equation 
in the symmetric nuclear matter is 

1 = C FU(u,k). (1) 

The factor Cq = Nq 1 , where Nq = 2ppm/Tr 2 , is the density of states for two sorts of the 
nucleons on the Fermi surface; F is the any of the Landau-Migdal dimensionless parameters of 
the quasiparticle-quasihole interaction F [9]: 



F(a x , n; a 2 , f 2 ) = C (F + F'{t x t 2 ) + G{a 1 a 2 ) + G'(tit 2 ) {B X B 2 )) 



(2) 



where a, f are the Pauli matrices in the spin and isospin spaces. The polarization operator 
U.{oj, k) is the integral over a nucleon particle-hole loop: 



U(u,k) = 4 



d 3 p 
(2^)3 



8(\p + k\- Pf)0(pf ~ p) 0(p - Pf)0{pf ~\p + k\) 



, p1 _ (p+fc) 2 , i 



2 m 



2m 



(p+fe) 2 . ■ 



(3) 

The factor 4 comes from the summation over the spin and isospin of the loop. The integral is 
taken by the traditional way [3]. The function k) can be represented in the form 



$(u>, fc) = fc) + <f)(-uJ, k). 



(4) 



The function <f>(u, k) is the Migdal's function [S] 

, . , . 1 m 3 \a 2 — b 2 
^ k)= 4^ 



In 



a + 6" 
a — b 



ab 



(5) 



where a = a; — (k 2 /2m), b = kpp/m^. 

Recall that the functions 4>(co, k) and 4>(—uj, k) (jSJ) have the following cuts on the complex 
w-plane. The cut of <f>(w, k) is placed at 



kpF k 2 
m 2m 



The cut of 4>(—oj, k) is placed at 



kpi 



2m 



< to < 



kp F k 2 
m 2m 



(6) 



< u < 



kpi 



m 



2m 



(7) 



1 This presentation of cj)(u!,k) differs from the one which was used in [17], but the results for the branches 
of solutions do not depend on the chosen form of <f)(ui,k) in (|4|5[) . In [17] we used the form which is more 
convenient in the study of the relation of the long wavelength instability to the pion condensation. 
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These overlapping cuts of the retarded polarization operator H(u, k) are shown in Fig.l. The 
polarization operator H(u, k) has the nonzero imaginary part, Im k)), for the frequencies 

u in the intervals ([6]), (171) . On the one hand, Im (II (u; , k)) appears as an imaginary part of 
the logarithm on the cut in ($5$). On the other hand, the same Im(n(a;, k)) can be interpreted 
as a result of the excitation of the real p/i-pairs for the given u, k, pp, when the function 
—i-n5{u) + t§^ — ^ ) contributes to integral ()3]). 

The physical sheet is defined by the condition that in the long wavelength limit k — > 
and v — > const, we obtain the physical sheet of the kinetic theory of Landau. As a result the 
imaginary part of the logarithm in Eq.(ll) in the complex w-plane is equal to (— iir) on the 
upper edge of the cut Eq. ([?)]) and (+iir) on the lower edge [2]. When the wave vector increases, 
the cuts overlapping becomes shorter and disappears at k = 2pp. At k > 2pp there are two 
nonoverlapping cuts in U(u),k). 

For not too large wave vectors k the collective excitations have the frequencies which are 
larger than the p/i-pair ones. It means that u s (k) lies on the real axis on the complex tu-plane 
more to the right than the cut fl6]). With increasing k, the cut overlaps with the solutions 0J s (k). 
The overlapping takes place at a specific wave vector k t h when the value of oj s {k th ) touches the 

k k 2 

cut Eq.([6]) at cu = th ^ F + Substituting uj s {k t h) in Eq.([T)) as solution we obtain 

2 + F ( x 
— = (x + 1 In 



x + 1 



where x = For F = 2 we obtain k th = 0.51 p (p =268 MeV). 

How can we continue our solution at k > fc^? At these wave vectors the dispersion equation 
(0Q) acquires the large imaginary part and we expect that the solutions are the imaginary ones. 
We are interested in the solutions with the negative imaginary part lmu s (k). There are two 
possibilities, the first: the branch u s {k) stays on the physical sheet and goes under the lower 
edge of cut and the second: u s {k) goes through the cut to the unphysical sheet. 

The logarithm In ( ^ m ^f PF ^if 2 ) = ln(Z) has the value ln(|Z|) — in on the upper edge of 



Lum—kpp + ^k 2 

the cut and the value ln(|Z|) + irr on the lower edge. Thus, on the physical sheet the 
discontinuity of logarithm over the cut is equal to (+2ni). Under the cut there is an unphysical 
sheet neighboring with the physical one on the Riemann surface of ln(Z). The magnitude of 
ln(Z) on the neighboring sheet differs by (— 2m). Thus, iri is canceled and we have a continuous 
change of logarithm, when go under the cut (following the arrow in Fig. 1). Then one may 
conclude that II (u;, k) in Eq. (CQ) changes continuously as well. 

To determine what of the possibilities realizes, we rewrite the Eq. ([!]) in the following way 



11 \ ■ 

— <f)(—uj, k) x 



C F(-A) 
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2 



1 -1 



In ( a -±4) . (9) 



We look for the solutions of Eq. (JTJ) and (Q with the negative imaginary part: u s = u r + kOi 
and LOi < 0. Let calculate the left (LHS) and right (RHS) sides of Eq.® at a definite k and Ui 
changing u r . In Fig. 2 the real and the imaginary parts of LHS is shown by the solid curves 



4 



(numbers '1' and '2'). The dashed curves are presented for the real and imaginary parts of 
RHS (curves '3' and '4'), the dotted curve is Im(RHS) on the nearest unphysical sheet (curve 
'5'). We see that the imaginary part of LHS (curve '2') is negative. The Im(RHS) (curve '4') 
is positive on the physical sheet under the cut. It becomes negative on the unphysical sheet 
(curve '5'). The solution is the crossing of the real parts of LHS (curve '1') and RHS (curve '3') 
and their imaginary parts (curve '2' and curve '4' or '5') at the same value of u r . In Fig. 2 we 
see the crossing of curves '2' and '5' therefore the logarithm must be taken on the unphysical 
sheet. The detailed discussion of Fig. 2 is done in Appendix. 

It is of importance that moving on the sheets of the Riemann surface of the one of the 
logarithms, we remain on the physical sheet of the other logarithm. The branch u s (k) goes 
under the cut fl6]) (Fig. 1) and not under the symmetric cut ([7]). We add a term (— 2ni) to the 
logarithm ln(Z) but not to another logarithm. To calculate u s (k) on the unphysical sheet, we 
solve Eq. ([T]) with the shifted logarithm. If we are interested in the symmetric branch oo' s (k) 
[Re(u' s (k)) = — Re (u s (k)) and Im (u' s (k))= Im (ou s (k))], then at k > k t h it goes through the 
symmetric cut ([7]). 

The branch u s (k) is presented in Fig. 3 by the solid curves. We see that oj s (k) is real at 
the wave vectors k < k t h and the imaginary part of u s (k) (shown in the bottom) appears at 
k > k t h- The dashed curves denote the cuts. The space between the x-axis and the lower 
dashed curve is the region where U(u,k) gets the imaginary part due to the overlapping cuts 
&, 0- The space between two dashed curves is the region where n(co>, k) gets the imaginary 
part due to the cut (Q. 

In this paper a simple model with the constant effective quasiparticle-quasihole interaction 
(T5]) is used. In Fig. 3 the solutions of Eq. (TjQ) are presented for the symmetric nuclear matter 
at the equilibrium density, corresponding to the Fermi momentum p j p=po=268 MeV. We take 
the quasiparticle mass m = 0.8mo (mo = 940 MeV) [9], [12] . The parameter F is a constant and 
F—2. This value of F corresponds to the isospin F' and spin-isospin G' parameters in the the 
Landau- Migdal effective interaction (recalculated to the value of C of this paper). 

2.1 Long wavelength limit 

In the long wavelength limit the dispersion equation ([T]) with Eqs. (T4|), ([3]) turns to the well- 
known Landau dispersion equation. For F > we can find k t h in ([8]). For k < k t h the solutions 
are real and close to the solutions of Landau dispersion equation. 

3 Dispersion equation with the real part of polarization 
operator 

It is instructive to compare u s {k) with the solutions of the widely used dispersion equation 
including the real part of the polarization operator only (7J El EE] : 

l = C FRe(U(uj,k)). (10) 
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The dotted curve in Fig. 3, having a "thumb-like" form, is the solution of this equation, u sr (k). 

We see that the upper part of the dotted curve and the solid one are identical when the 
polarization operator is real. But at k > k t h there is a difference between them. The mode 
u) s (k) becomes complex- valued. The real part Re (u s (k)) continues to increase. The imaginary 
part lm(u s (k)) is negative and smoothly decreases. The solutions of Eq. ljTOl) demonstrate 
another dependence: u) sr {k) decreases at k > k th and disappears at k = k max = 0.74 p . The 
attempts to find additional solution of the dispersion equation ([T]) (i.e., with the full polarization 
operator) corresponding to the lower part of the dotted curve u sr (k) (Fig. 3) on the physical 
or neighboring unphysical sheets failed. 

When the damping is small, we can use the following equation for the imaginary part of the 
zero-sound mode [3]: 



Here u sr (k) are solutions of Eq. ffTUj) . which are placed on the upper part of the "thumb-like" 
form in Fig. 3. In Fig. 3, u S i(k) is shown by the dotted curve. We see that u S i(k) and Im (u s (k)) 
coincide at a very small damping only. 

In some papers in which the pion dispersion equations are considered [Hi HE], the equations 
with the real part of the polarization operator (similar to Eq.f JTO]) ) are used. The solutions 
have a complicated form with a lot of the "thumb-like" structures. The method suggested in 
this paper, permits us to obtain several one- valued smooth complex functions 0J n {k) instead of 
the many-valued "thumb-like" structures. Moreover the branches of the solutions u) n (k) can be 
related to the solutions on the physical sheet by the uninterrupted continuation on k or p F or 
F |17j . Thus we know the physical nature of the every branch of solutions. Since the imaginary 
part of the solutions appears due to the passing to the unphysical sheet through the cut, the 
nature of the cut determines what kind of the width of excitations we calculate. This statement 
is valid, for example, for the zero-sound excitations in the nuclear matter in the spin-isospin 
channel when both the particle-hole and the isobar-hole cuts exist. 

4 Giant resonances 

In this Section the giant dipole resonances are considered. In the simple version of the Migdal 
particle-hole effective interactions the electrical dipole collective states are excited by the isovec- 
tor interaction constant F' (J2J). The coefficient Cq is defined in (JTJ) . The values of F' are changed 
in the region from 0.8 to 1.9 Q2] at C = 7r 2 /(2po™o) = i50 MeV fm 3 . 

First of all we obtain the solution of Eq.flTJ with F = F' . The Steinwedel- Jensen model 
[TU [15] determines the wave vector kA which corresponds to the definite excitations in the 
nucleus. We follow to |15j . where = where R = tqA 1 ^ 3 , Vq — 1.2 fm. Using this relation 
of k to A 1 / 3 we get u s (A). We see that k is larger for the lighter nuclei (this is the consequence 
of the boundary conditions). 

In the next step we compare u s (A) with the formula for the GDR energies in the different 




(11) 
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nuclei [16] : 



E GDR (A) = 31.2A" 1/3 + 20.6A- 1/6 MeV. 



(12) 



Before doing the comparison we investigate the dependence of the branch u s (k) on the 
nuclear matter parameters: pp, F, m which are included in Eq. ([1]) . In Fig. 4 we consider the 
dependence of u s (k) on the constants of the effective quasiparticle interactions F (Fig. 4a), on 
the value of Fermi momentum pp (in Fig. 4b) and on the effective mass m (Fig. 4c). 

In Fig. 4a we present the branches u s (k) for that values of F which are obtained from 
the experiments [12]. In [3] it is demonstrated that the frequency of zero-sound excitations 
increases with increasing of F (Fig. 16.1 [3]) at a fixed k. We obtain the same behavior if to fix 
rather small k and take for F such values that u s (k) is real. But in the region of overlapping 
of the zero-sound and particle-hole modes, when Imu s (^) is nonzero, the behavior is changed, 
this is shown in Fig. 4a. The growth of F results in decreasing of both R,eu s (k) and lmu s (k). 
The sensitivity of the real and imaginary parts of u s (k) to the value of F is different: Keu s (k) 
is almost equal at F — 1 and F = 2 but lmu s (k) is changed considerably. 

In Fig. 4b we present the branches u s (k) for some values of pp. The growth of pp from 
200 MeV to 300 MeV results in the increasing of R,eu s (k) and decreasing of lmu s (k). In Fig. 4c 
we see that the growth of the effective quasiparticle mass from m = 600 MeV torn = 940 MeV 
give the decreasing of both Recu s (k) and lmuj s (k). 

Thus, we obtain a different dependence of the zero-sound excitations on the nuclear matter 
parameters F', pp and m. Changing these parameters we try to approximate the phenomeno- 
logical equation fTT2l by u) s {kA)- 

In Fig. 5 the comparison of the experimental and calculated energies of the giant dipole 
resonances in the nuclei is presented. It is difficult to expect that such a complex phenomenon 
as giant resonances can be described quantitatively by the simple model used in this paper. 
Nevertheless, we expect that the main features are reproduced. In Fig. 5 we present the curve 
(TT2"j) by the solid line. The Reo; s (kA) correspond to the energies of GDR and are shown at the 
positive u>. The Imo; s (kA) which correspond to the half-width of the GDR are negative and are 
shown at the negative u. The dashed-dot curves show the results for the same nuclear matter 
parameters as in Fig. 3 (pp = 268 MeV, F=2.0, m = 0.8m ). We see that the calculated 
energies of resonances are too large and the width too small. 

Now we try to find that parameters of the nuclear matter which permit us to describe, at 
least generally, the GDR in the nuclei. To approximate the phenomenological curve f)12p by 
our results we use two sets of the nuclear matter parameters. One set is the more suitable for 
the light nuclei and the second is better for the heavy ones. The branch u s (kA) obtained with 
Pp = 200 MeV, F—1.2, m = mo is shown by the dashed line. By the dotted line we draw the 
u) s (kA) obtained with the matter parameters pf = 200 MeV, F—1.6, m = 0.8mo- 

For the nucleus 12 C we have from Eq. fTl2|) Egdr = 27.2 MeV, the experimental escape 
width is T w 3 MeV. On the dashed curve we obtain Egdr = 30 MeV and T = 6.8 MeV. On 
the dotted curve: Egdr — 37 MeV and T = 7.8 MeV. For the nucleus 208 Pb we have from 
Eq. (fT2"l) Egdr — 13.7 MeV, the experimental escape width is T w 0.5 — 2.0 MeV. On the dashed 
curve we obtain Egdr = 10 MeV and T w 0.5 MeV. On the dotted curve E GDR = 12.7 MeV 
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and T = 0.5 MeV. We obtain the better approximation to the curve Eq. (jT2l) when start from 
the nuclear matter with rather low density: pf=200 MeV, p/po ~ 0.42. This points to the 
important contribution of the nuclear surface to the excitation of the resonances. 

In Fig. 5 the decreasing of the GDR energies and escape width with the growth of the 
atomic number A is demonstrated. In Fig. 3 we see the expected figure when the Reu; s (k) 
and (— Ima; s (k)) increase with the growth of the wave vector k. The Steinwedel- Jensen model 
[T4"l [15] shows that the excitations with the same quantum numbers have the larger wave vectors 
k in the lighter nuclei. Thus using the simple model of the zero-sound in the nuclear matter 
we obtain the qualitative description of the experimental fact that the energies and the escape 
width of the giant resonances in the light nuclei are larger than in the heavy ones. 

5 Discussion 

In this paper we demonstrate the method of calculation of the zero-sound frequencies (real or 
complex) in the region of overlapping of zero-sound and particle-hole modes. It is shown that 
the damped zero-sound solutions are placed on the unphysical sheet which is neighboring with 
the physical one on the Riemann surface of the logarithm in Eq. ([5]). The logarithmic cut fl6]) 
corresponds to the excitation of p/i-pairs, and we can conclude that Im (uj s (k)) is the half-width 
of the zero-sound mode decay into p/i-pairs. 

It is shown in the Fig. 3, that neglecting of the imaginary part of the polarization operator 
in the dispersion equation, we obtain the branch of solutions having a "thumb-like" form with 
the upper point of the real part and the quick decrease of the imaginary part at k = k max . 
On the other hand, considering the full polarization operator at a proper unphysical sheet, we 
have an continuous behavior of the zero-sound mode oj s (k) without the singularities. 

The zero sound branches u s (k) are calculated in the nuclear matter with the different values 
of the Fermi momentum pp, the nucleon effective mass m and the effective residual interaction 
F. The giant resonances can be considered as the zero sound excitations in the nuclei. Using the 
simple model of the quasiparticle interaction ([2]) we obtain at least the qualitative description 
of the experimental fact that the energies and the escape width of the giant resonances in the 
light nuclei are larger than in the heavy ones. 

The zero sound mode uj s (k) calculated in the nuclear matter in equilibrium (pp = 268 MeV, 
F=2.0, m = 0.8m , C*o=150 MeVfm 3 in Eq.([T])) corresponds to the very high energies and low 
escape widths of GDR (dot-dashed curve in Fig. 5). The GDR in the nuclei correspond rather 
to the zero sound in the nuclear matter with lower density (the dashed and dotted curves in 
Fig. 5). This fact is probably related to the importance of the nuclear surface in forming of the 
resonances. 
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6 Appendix 



In this Appendix we try to show that the solutions of the zero-sound dispersion equation are 
placed on the nearest unphysical sheet when the overlapping of the zero-sound and the ph 
excitations take place. 

In Fig. 2 the left and right sides of Eq. (Q are shown as functions on the frequency. The 
calculation is made for the fixed wave vector k=0.Q5pQ. As it is shown in Fig. 3, the solution of 
Eq.QTJ) (and of ©, identical to it) at k=0.65p is equal to u s (k = 0.65p ) = (0.304, -0.003)p - 
We choose fc=0.65po which is close to k t ^ in Eq. (|SJ), and we see in Fig. 3 that the solution of 
§B) is close to the solution of the approximate equations (ITU]) . (TTTT) . 

In Fig. 2 the imaginary part of u is taken to be equal to the imaginary part of the solution 
Ui =-0.003po- The solid curve '1' ('2') is the real (imaginary) part of LHS of Eq. (J5J). We see 
that for the negative Wj the imaginary part of LHS is negative also. The dashed curve '3' is the 
real part of logarithm in the RHS of (Q and the curve '4' is its imaginary part on the physical 
sheet. As it was discussed above the imaginary part of logarithm on the physical sheet on the 
lower edge of the cut is positive and close to +iir. The dashed curve '5' is the imaginary part 
of logarithm on the nearest unphysical sheet of the Riemann surface, where the imaginary part 
is differ by —2m and is negative. In Fig. 2 a letter 'A' ('B') denotes the crossing of the real 
(imaginary) parts of LHS and RHS of (jUJ) at u r = 0.304p , that corresponds to the solution of 
([T]). In Fig. 2 we see that the imaginary parts of LHS crosses that branch of logarithm in RHS 
which is place on the unphysical sheet. 
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7 Figure captions 

FIG. 1. Complex plane of the frequency u. The cut fl6]) is shown from 1 to 1', the cut (J7J) is 
shown from 2 to 2'. 

FIG. 2. Comparison of LHS and RHS of Eq. ©. Solid curves: '1' presents Re (LHS) and '2' 
shows Im(LHS). Dashed curves: '3' is Re (RHS), '4' is Im (RHS) on the physical sheet, '5' is 
Im (RHS) on the unphysical sheet. The figure is made at k — 0.65po and Ui = — 0.003po- A letter 
'A' ('B') denotes the crossing of the real (imaginary) parts of LHS and RHS at u r = 0.304p - 
p =0.268 MeV. 

FIG. 3. Zero-sound branch of solutions u s (k) of Eq. (JTJ (solid curves). The u sr (k) are the 
solutions of Eq. f JTOj) and u S i{k) are solutions of ffTT]) (dotted curves). Dashed lines are the 
boundary of the cuts (Q, ([7]); k t h denotes the wave vector when the overlapping of u s {k) and 
quasiparticle continuum starts. 

FIG. 4. Zero-sound branch of solutions u s (k) for the different values of the effective constant 
F (figure a), the Fermi momentum pp (figure b) and the effective mass m (figure c). Re (u s ) 
are drawn at the positive u, Im (u s ) are drawn at the negative u. The blobs mark k t h/po f° r 
the curves 1,2,3. Fig. a: F—0.2, 1.0, 2.0 (numbers 1,2,3, correspondingly). Fig. b: j9p=200, 
268, 300 MeV (numbers 1,2,3, correspondingly). Fig. c: m=500, 750, 940 MeV (numbers 1,2,3, 
correspondingly) . 

FIG. 5. Comparison of the experimental GDR energies with the calculated ones (A is the 
atomic number). At the positive co we compare Egdr (fT2l with Reco> s (A). At the negative u 
the Imo; s (A) are presented. The solid curve shows the approximation (fl2"j) . Dashed curves stand 
for u 8 (A) calculated for p F =200 MeV, m = m and F=0.6. Dotted curves are calculated with 
Pi?=200 MeV, m = 0.8m and F=0.5. Dash-dotted curves are for p F =268 MeV, m = 0.8m 
and F=2. 
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